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A METHOD FOR THE DESIGN OF POROUS-WALL WIND TUNNELS 
By George M. Stokes 


SUMMARY 


This report contains a description of a general method for the 
design of a porous wall (slotted and perforated) wind tunnel which can 
be applied to the design of low Mach number supersonic wind tunnels . 

This design method is concerned with the development of a uniform super- 
sonic flow in the wind tunnel and the selection of the proper test-section 
open ratio necessary to minimize shock- wave-ref lection interference. The 
main feature of the method is its ability to indicate how design require- 
ments and wall porosity characteristics can be combined to achieve a good 
flow for desired conditions. A limited amount of experimental data has 
been obtained to show that the method can successfully produce a suitable 
flow field for model testing. More research is necessary, however, before 
the practical significance of the method with regard to wave-reflection 
cancellation can be determined. 

The experimental data obtained from a perforated-wall wind tunnel 
designed by this method showed that the test-section Mach number varia- 
tion along the center line was ±0.007 for the design Mach number, 

M = 1.278. 

A system for calculating the flow in a porous-wall wind tunnel has 
been developed. This technique is expected to be useful also in deter- 
mining the flow in the test section at off -design supersonic Mach num- 
bers. When this system was used to calculate the flow in a 3- by 3-inch 
wind tunnel, it was found that the construction errors in the wall open 
ratio were sufficiently small to cause no adverse flow variations. 

A discussion is included to explain why the open ratio as given by 
the design method used in selecting the open ratio of the wall does not 
prevent a disturbance from being generated at the point on the wall where 
the model shock wave intersects. 


INTRODUCTION 


In the interest of developing methods for the design of transonic 
tunnels having slotted and perforated walls, several different studies 
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have been made (e.g., refs. 1 and 2). Reference 1, which contains a 
simplified analysis of the tunnel-wall characteristics required to gen- 
erate shock-free supersonic flow in minimum distance, indicates that a 
constant -porosity wall selected on the basis of calibrations with the 
flow normal to the wall will not generate the desired flow. It does sug- 
gest theoretically, however, that a suitable flow distribution could be 
developed by the use of porous walls if the walls would remove flow at a 
specified rate at each tunnel station. The study of transonic -wind- 
tunnel design presented in reference 2 goes a step further in the devel- 
opment of a flow generating method by considering slot discharge coeffi- 
cients for relating the wall outflow to the slot width. Although in 
general the principles proposed in these previous design investigations 
are sound, the work is incomplete to the extent that it does not indicate 
the influence or treatment of the boundary layer, nor does it show how 
to determine the required wall-open-ratio distribution to satisfy desired 
flow conditions. 

Since it is important in the cases of shock-wave-reflection studies 
(refs. 3 and 4 ) to have a design method available whereby certain desired 
boundary conditions may be stipulated, further study has been devoted to 
the development of a porous-wall transonic-design method. The problem 
was investigated from the points of view of keeping the flow-generation 
region short, thus obtaining minimum power loss, and of controlling the 
thickness of the boundary layer at the end of the flow-generation region. 
The results of this study yielded a successful design method. This 
report contains a detailed description of the developed method. 

The general design method presented herein covers a method or sys- 
tem for finding the open ratio required to minimize the shock-wave reflec 
tion for constant boundary-layer growth conditions and a method for 
obtaining within limits any chosen tunnel-flow distribution for a two- 
dimensional wind tunnel. To show the effectiveness of the design method, 
experimental results are presented for a perforated-wall wind tunnel 
designed by the method. 


SYMBOLS AND DEFINITIONS 

a speed of sound in air, ft /sec 

Cf effective wall friction coefficient (ref. 3) 

H total pressure, lb/sq ft 

Hj total pressure in tunnel main stream over porous wall, 

H]_ - Shock losses, lb/sq ft 
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h 

height of tunnel, in. 

o 

M 

Mach number 


“nt 

total mass-flow rate per unit 
slugs/ft^-sec 

area through porous material. 

m nt ' 

total mass-flow rate per unit 
for 6=0, slugs/ft^-sec 

area through porous material 

“nbi 

mass-flow rate per unit area through porous material required 
to compensate for boundary- layer growth rate, slugs/ft^-sec 

’“ex 

component of stream mass-flow rate normal to center line pro- 
viding for supersonic expansion, slugs/ft^-sec 


P static pressure in stream, lb/sq ft 

p. static pressure in plenum chamber (ref. 5) > lb/sq ft 

J 

AP static pressure drop across porous material, Pg - Pj , 

lb/sq ft 

T temperature 

V velocity, ft /sec 

V nt component of velocity normal to porous material, i%t/ p s> 

ft /sec 

Vj jet velocity, velocity corresponding to pressure ratio Pj/Hj 

across porous material (ref. 5)> ft/sec 

r open ratio, ratio of open area of porous material to total 

area of material 


w width of wall, in. 

6 stream deflection angle caused by model, or stream angle with 

respect to tunnel center line (5 positive for streamlines 
diverging from center line), deg 

6* boundary- layer displacement thickness, in. 


dS*/dx 


rate of change of boundary- layer displacement thickness 
over wall 
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d0/dx rate of change of boundary-layer momentum thickness over 

porous wall 

p density of gas before entering porous material, slugs/cu ft 

Pj density of gas in jet (ref. 5), slugs/cu ft 

p stagnation density behind oblique shock 

v angle through which a supersonic stream is turned to expand 

from M = 1 to M > 1 (ref. 6), deg 

v a stream expansion angle resulting from flow expansion of the 

same family or same wall, deg 

stream expansion angle resulting from expansions of family 
opposite to that of v a , deg 

a discharge coefficient for parallel flow, ^t / p jVj r (see 

ref. 5) 

y porous -wall convergence angle (angle positive with walls 

converging downstream toward tunnel center line), deg 

7 1 ratio of specific heat of gas at constant pressure to spe- 

cific heat at constant volume 

Subscripts: 

0 stagnation conditions 

1 . free-stream conditions in test section 

2 conditions behind oblique shock caused by model 

ah ahead 

j nonviscous jet conditions with plenum-chamber static pressure 

(ref. 5). and stream total pressure 

bZ boundary layer 

p porous wall 


max 


maximum 
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pt point 

s stream condition 

sw solid wall 

The word "porous" is considered throughout the paper to apply to 
both perforated and slotted materials. 

The words "plenum chamber" are used herein to denote the low-pressure 
chamber on the opposite side of the porous wall from the test section. 


GENERAL DESIGN METHOD 


The general design method has been developed to. apply to two- 
dimensional wind tunnels having porous top and bottom walls and solid 
side walls. (See fig. 1.) As shown in this figure, the porous walls 
are vented into a plenum chamber which is attached to a suction supply. 
The upstream portion of the porous material is used for generating the 
required supersonic flow by means of expansion associated with outflow 
into the plenum chamber. The remaining porous material is used to form 
the walls of the test section. 

In the design of transonic tunnels, it is necessary to consider the 
wall interference effects resulting from reflecting disturbances as well 
as the uniform requirements of the test section. Because of this the 
general design method is presented in two parts; one part shows how to 
determine the open ratio needed in the test section to minimize the shock 
reflections, and the other shows how to determine a wall -open- ratio dis- 
tribution in the flow-generation region which will give a uniform flow 
in the test section. Since an important consideration is given to the 
boundary-layer growth by the design method, a brief description of the 
boundary-layer behavior with respect to porous materials will be pre- 
sented first . 


Estimation of Boundary-Layer Growth 

Only a little published information appears to be available (refs. 3 
and 7) to show how the boundary layer grows over porous materials or what 
the friction coefficients for such materials might be. Recently, however, 
a limited number of tests were made in the Langley. 3- by 3-inch transonic 
flow apparatus to determine the boundary- layer growth rates over a 

4l-percent-open porous specimen. The test specimen was a 2^- inch-square 

perforated sheet 0.060 inch thick, having 0.0407-inch-diameter holes. 
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Pressure measurements were made to determine the thickness of the bound- 
ary layer at two points^ one point was l/8 inch upstream of the specimen 
and the other was l/8 inch downstream of the specimen. From these data, 
the boundary-layer growth rate d&*/dx p was found by dividing the dif- 
ference between the downstream and upstream boundary-layer-displacement 
thickness by 2^ inches. The results of the tests are shown in figures 2, 

3, and 4. 

The data of these three figures were taken with the upstream boundary- 
layer displacement thickness, 5*, equal to approximately 0.002 inch, and 
with the stagnation constants H 0 and T 0 equal to 30 inches of mercury 
and’ 200° F, respectively. Figure 2 shows how the boundary-layer growth 
rate d6*|dXp changes with the outflow velocity ratio V nt jV s for several 

values of the free-stream Mach number. These data are presented in this 
form to expedite their use with the design method. Curves describing the 
boundary- layer growth relationship should be obtained for a range of open 
ratios to define completely the growth rate with respect to the open ratio. 
Thus, the data presented in figure 2 are inadequate for determining accu- 
rately the boundary- layer growth rates over similar specimens at differ- 
ent open ratios, and for other upstream boundary-layer thicknesses. Never- 
theless, if some approximate calculation is desired, it is believed that 
these data may be used for estimating the boundary-layer growth rates for 
open ratios within 10 to 20 percent of the measured value of open ratio. 

One method used for making these estimates is to assume (l) that a linear 
relationship exists between the open ratio and the rate of growth of the 
boundary layer at zero outflow and (2) that the rate -of -growth curves for 
other open ratios would have approximately the same shape as a function 
of outflow velocity ratio V nt /V s as the curves for the 4l -percent-open 

material. By use of these assumptions, the data of figure 2, and a solid- 
wall boundary-layer growth rate of 0.002 (ref. 3), the boundary- layer 
growth rate over a 20-percent-open perforated plate would be estimated to 
be 0.004 for a zero outflow velocity condition at a Mach number of 1.01. 

To estimate the total boundary-layer growth rate in wind tunnels having 
both solid and porous walls, the boundary-layer growth of the solid wall 
must be added to that of the porous walls. The following expression 
shows how these two boundary- layer growth conditions enter the porous- 
wall outflow equation: 


m nt 


Ps^nt 

p s V s sin(& + 7) + P S V S cos(5 + 7) + P s ^s 


h sw d5* 
Wp dx sw 


cos 7 


( 1 ) 
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This equation was obtained by assuming the flow conditions pictured 
in the sketch below. 

-Velocity component parallel to porous 
wall, equal to V s cos (5 + 7) 

Line parallel to tunnel center line 

Velocity component normal to porous 
wall, equal to V s sin(6 + 7) 

Direction assumed for flow 
between porous wall and 5* 



arc tan 


Effective free-stream 
boundary 


d5* 

dxr, 


Porous wall 


Line parallel to 
center line 


The term p s V s sin(S + 7) of equation (l) denotes the mass rate 
of outflow through the porous wall resulting from the flow angle 8 


and the wall convergence 
h 


The terns 


P S V S cos (§ + 7) a nd 
dx p 


P S V S - cos 7 of equation (l) represent the additional mass rate 


dx 


"p "sv 

of outflow taken through the porous wall to account for the boundary- 
layer growth over the porous walls and side walls, respectively. In the 
development of equation (l), it was assumed that the flow crossing the 
effective free-stream boundary over the porous wall turned in a direc- 
tion normal to the wall and passed through. This consideration is equiv- 
alent to assuming that a void exists between the effective free-stream 
boundary and the wall. The last term of equation (l), representing the 
mass rate of flow through the porous wall needed to compensate for the 
increase of boundary-layer displacement thickness over the side wall is 
an approximation, since in general the flow angle on the side walls is 
different from zero and variable over the wall height h sw . For con- 

d 5 * 


venience, the sum of the boundary- layer growth terms p_V 


s v s 


dx r 


cos (5 + 7) 


and p s V s -S* 


hsw d&* 


w^ 


”p ^ x sw 

of uniform flow in the wind tunnel (i.e. 


cos 7 will be denoted as 


For the condition 


m nbl- 

8=0), equation (l) becomes 


“nt' = P s V s( tan 7 + g- + 


d&* , hsw d&* 




dx 


sw 


cos 7 


( 2 ) 
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It should be noted that the solution of equation (l) for m n t can- 
not be obtained directly. This results because the boundary- layer growth 
rate over the porous portion of a wind tunnel is a function of the total 
outflow m n t and also of V nt ^V s . if all the terms of equation (l) are 
known except the term containing d5*/dxp, a solution of the equation 


can be found by an iteration process. This may be accomplished by first 

assuming = 0 and then finding the value of the V nt /v s given by 

dx p 

the other terms of equation (l) . This value of V^/Vs is used to deter 

mine (from curves for the particular material similar to those of fig. 2 ) 
an approximate value of d5*/dXp to be tried in the next calculation for 

Vnt/Vg. This iteration process is repeated until the value of V n t/V s 
resulting from equation (l) is equal to the value used to read d 6 */dx p . 
The value of m n t obtained by use of this converged value of d5*/dx p 
is the required result. 


When using the data presented in figure 2, it should be remembered 
that the boundary- layer growth over a perforated material depends not 
only upon the velocity ratio V n ^|V s and Mach number, but also upon 
other flow variables. This can be' shown if the following form of the 
momentum equation 



Ps v nt _ £f 
Ps V s " 2 


is rewritten in terms of boundary-layer displacement thickness. This 
gives 


d&* = &* dHi _ 5* dM 

(2 - M 2 ) + H' 

+ H' 

PsVnt . Cf 

dxp H ' dx M dx 

7 f - 1 

1 + 21 — M 

P S V S 2 


where H' is the ratio of the boundary-layer displacement thickness to 
the momentum thickness and 7 ' is the ratio of the gas specific heat at 
constant pressure to the specific heat at constant volume. Here, it is 
seen that the boundary- layer growth is a function of the displacement 
thickness, the form factor H* , the stream Mach number and velocity, the 
porous wall outflow velocity, the wall friction, dM/dx, and dH'/dx. 

The most significant point to be made from the data of figure 2 is 
that the boundary-layer growth rate drops rapidly as the porous -wall 
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outflow velocity is increased and that the stream Mach number has an 
appreciable effect on the growth rate. 

Figure 3 shows how the effective friction coefficient (for uniform 

flow) + s n (eq. (3) of ref. 3) varies with outflow velocity 

2 dx P Ps v s 

and Mach number. The data are presented primarily to show the approxi- 
mate magnitudes of the friction coefficients for perforated material. 
Curves of this type are not generally needed when data are available 
which show the growth of the boundary layer. In addition to the boundary- 
layer data taken for the 4l-percent-open specimen, the porosity charac- 
teristics of the material were also measured. The curves of figure 4 
show the relationship between discharge coefficient a and velocity 
ratio V s /Vj for several Mach numbers. The dashed portions of the 
curves for the Mach numbers of 1.17 and 1.27 have been extrapolated 
from the basic data shown. This extrapolation was necessary in order 
that the range of the data would be sufficiently large for use in the 
design of the perforated-wall wind tunnel described later in the paper. 


Method for Selecting the Test -Section Wall-Open Ratio 

It is well known that shock-wave reflections from wind-tunnel walls 
in the transonic speed range are quite troublesome. Presently, much 
research is being conducted toward the development of porous walls in an 
attempt to reduce the effects of wave reflection. The method presented 
in this section deals with the selection of the porous-wall open ratio 
needed in the wind-tunnel test section to reduce or minimize the shock- 
wave reflecting ability of the tunnel walls. 

The method is based on the reasoning that the porous -wall boundary 
of the wind tunnel will not reflect a model shock- wave disturbance in 
the stream if the normal component of the stream flow is absorbed by the 
wall; that is, the flow approaching the wall must be completely removed 
at its incoming rate. Thus, the problem is one of determining the wall- 
open ratio necessary to fulfill this condition. 

An equation developed in reference 5 is 


r = ■ m nt 

‘ " p jV 


(3) 


This equation shows how the wall-open ratio is related to the wall out- 
flow, the theoretical jet conditions pj and Vj, and the wall porosity 

characteristic a. Equation ( 3 ) is used in combination with equation (l) 
to calculate a value of open ratio for the walls of the test section. 
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Equation (l) is rewritten for use in this section as 


m nt2 = P2 V 2 


sin 


(So + 7) + cosfsp + 7) + 4^— cos 7 

\ / dx B . V <= ) Wp dx sw 


00 


because the conditions (subscript 2) downstream of an obliqud shock are 
being considered. Similarly, the quantities pj, Vj, and m nt in equa- 
tion (3) are defined as Pj2> v j2> an<i m nt2 when used to refer to the 
conditions downstream of the oblique shock. 

When solving equations (3) and (4) to find the test-section open 
ratio it must be assumed (l) that a uniform supersonic flow enters the 
test section, (2) that the porosity characteristics of the wall type 
intended for use have been defined similarly to those of figure 4, 

(3) that the side-wall boundary- layer growth rate is known, (4) that 
boundary-layer growth rate data similar to those of figure 2 are avail- 
able, and (5) that the wall convergence angle has been chosen. Equa- 
tion (4) should first be solved for m n t2 for several arbitrarily chosen 
values of 6 which correspond to values expected at the walls of the 
tunnel when tests are being conducted. 

The values of r which will result from insertion of these values 
of m n t2 in equation (3), will show that a slightly different value 
of r is obtained when 6 is changed, thus the value of r to be 
chosen, should be a value which corresponds to a value of 8 that best 
represents the range of deflection angles expected during wind-tunnel 
tests. The values of Pj 2 an( i V j2 ^ use( i in equation (3) must 

result from use of a value of the plenum-chamber pressure required to 
provide a uniform flow in the test section for the value of test -sect ion 
open ratio which is being determined. Since in equation (3) the quanti- 
ties pj 2 > r > and m nt2 (because of d8*/dxp^ are dependent var- 

iables, the solution of the equation must be obtained by use of a cut- 
and-try process . 

Because of the indirect method which must be used to obtain the 
test -section wall open ratio and because of the complexities involved 
in solving equations (l) and (3), a detailed step-by-step procedure is 
introduced here to expand upon the summary solution of the problem just 
stated and to indicate a practical solution. This procedure is as 
follows : 


1. Choose the desired test-section Mach number, M^. (The chosen 
Mach number is expected to be in the supersonic range M = 1.00 to 
M = 1.50.) Also stipulate the tunnel stagnation conditions H 0 i and 
T 0 i, and the convergence angle 7. 
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2. Determine from oblique shock theory (ref. 6 ), the stream condi- 

tions H 2 /H o1 , P2 /p o 1^ a 2 / a ol^ Po2/P 0 l^ and m 2 downstream of the 
shock wave for several values of b 2 from 0 to & 2 = b pma x * If is 
advisable to plot these quantities versus 5 2 . This will account for 
all values of strengths of the weak shock which may be encountered in 
supersonic flow up to • Expansion waves are not accounted for by 

this design method. 

3. For the stagnation conditions of step 1, and by the use of ref- 
erence 6 , compute H 2 , p 2 , V 2 , and p o2 for the range of S 2 f s chosen 

in step 2 . 

k. Neglect the boundary- layer terms and the convergence angle 7 
in equation (A) for the present and determine the resulting values of 

m nt 2 g iven 


m nt 2 = P 2 ^ 2 s ^ n ^2 


to obtain approximate values of the required test-section mass-outflow 
rates. Use the values of p 2 and V 2 from step 3 for the range of S 2 
in step 2 . 


5. Neglect the boundary- layer growth in the test section ahead of 
the shock and the wall convergence angle 7 so that P]_ = P j . Also 

take the total pressure Hj 2 = H 2 and compute Pj|Hj 2 , for the range 

of 5 2 *s in step 2 . 

6 . From the pressure ratios Pj^Hj 2 obtained in step 5* calculate 

by use of reference 6 the corresponding values of p^ an< ^ Vj 2 (assuming 


isentropic flow) using the relations 
where P oj2 = p q2 and a oJ2 = a ol = 


p j2 

Pj 2 = p oj 2 > 
h oj2 

7'RT ol . 


V j2 - Mj 2 


a j2 

a oj 2 


a oj2 


7 . The quantity a needed in the evaluation of equation (3) must be 
determined from curves showing the porosity characteristics of the wall 
type for use. If no such data are available, it will be necessary to 
determine a experimentally for a sample of material having a value of 
open ratio guessed to be in the range required. The data of reference 5, 
which show some curves of a plotted against V s/Vj> indicate that it is 

not* essential to make the calibrations for determining a for the exact 
open ratio of the material required in the final configuration. 
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8. By use of ni n t2 s ^ e P Pj2 ant ^ ^j2 s ^ e P and values 
of a for V2^Vj for the corresponding M2, calculate using equation (3), 
the approximate values of open ratio required for the range of 62' s chosen 
in step 2. 

9. Select from the open-ratio values obtained in step 8 a value of 
open ratio which corresponds to a flow deflection angle, &2, which best 
represents the range expected during wind-tunnel tests, and determine 
from experiments on a sample of material having this value of open ratio 
curves of d6*/dxp as a function of V n t/V s for several values of M 

(in the range of M2). Also determine curves of a similar to those of 
figure 4 for this material to check the values previously used. 

10. Start with this step to recalculate the values of open ratio 
required in the test section for various values of 62 with the tunnel 
boundary layer included. Estimate d8*/dx sw by use of the momentum 

equation in conjunction with reference 8, or by use of experimental 
measurements . 

11. Using the curves obtained in step 9 and the dS*/dx sw esti- 
mated in step 10, obtain m^' (by the iteration procedure previously 

discussed in conjunction with eqs . (l) and (2)) for the region of paral- 
lel flow in the test section upstream of the shock. Also obtain the 
values of m nt2 for the range of 82 in step 2, including the boundary- 

layer terms and the convergence angle 7 in equation (4). 

12. In order to find the plenum-chamber pressure with the boundary 
layer included in the calculation, proceed by the use of reference 6, 
and find Mjp, Pji^P Q ji^ and - a jl/ a ojl for a chosen range of pressure 

ratios slightly less than the pressure ratio value given by 

Pi/Hqi, and calculate values of 


-ntl = = or gi- P olMjl |ii- a 0l 


(where P 0 ji - P Q i and a 0 jp = a Ql ^ for this range of pressure ratios 
using the r selected in step 9. 


J Sji 


H„ 


13- Compute = - — n p 

chosen in step 12 and plot p, 

J 

of determined in step 12, 


(where = H ol ) for the values Pj/H^ 

as a function of the corresponding values 
Read a value of Pj from this 


curve 
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for a value m nt i equal to the value of m n t' determined in step 11. 
This gives the plenum-chamber static pressure required to provide uni- 
form flow in the empty test section. 


1^. Using the Pj determined in step 13 and Hj 2 = H-2> recalculate 
Pj/Hj2 and the corresponding values of pj 2 and Vj 2 , as in step 6. 

15. Make new calculations for r (eq. (3)) using the values m n -^ 2 
from step 11, with the corresponding values of pj 2 and Vj2 from 
step 14, and the curves a as in step 8 for the range of 5's considered. 


The newly calculated required open-ratio range resulting from the 

inclusion of the boundary- layer growth will in most cases indicate that 

a test-section open ratio different from the one estimated in step 9 (no 

boundary 'layer) is needed. Thus, several iterations (using steps 9 to 15) 

may be required before the final open-ratio selection made is equal to the 

value used for estimating the boundary- layer growth and finding P<. 

J 


The above method has been found to be effective in determining ana- 
lytically the required open ratio of a porous wall. It must be borne in 
mind, however, that as indicated, the method accounts for only one boundary- 
layer growth rate along the wall behind an oblique shock for a given 6 
and y; thus, if the growth rate changes appreciably with distance for 
different downstream stations through this region, a constant value of 
open ratio in the test section will not provide the necessary outflow. 


Flow Generation 

This section describes a method for designing a porous supersonic 
wind-tunnel nozzle capable of providing a uniform supersonic flow at a 
chosen Mach number in a given distance. It is assumed when using this 
method that a uniform sonic flow will enter a two-dimensional supersonic 
flow-generation section of the tunnel, and that the walls providing the 
flow expansion will be porous and will be surrounded by a plenum chamber 
in which the pressure can be controlled. The schematic drawing in fig- 
ure 1 pictures a representative tunnel configuration. 

To initiate the design of the nozzle, certain quantities must be 
chosen; these are the test-section Mach number, the tunnel stagnation 
pressure and temperature, the test-section open ratio r, the rate of 
supersonic flow expansion for v a family, the convergence angle, the 
physical dimensions of the tunnel, and the side-wall boundary-layer 
growth rate db*/dx sv . With this information and a requirement that the 
porous wall originate at the minimum area or point where M = 1.0, a two- 
dimensional supersonic characteristic net (ref. 9) must be constructed 
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for the chosen expansion rate, or center-line Mach number distribution, 
to obtain a theoretically desired flow for the flow-generation section of 
the tunnel. From this characteristics network, the flow quantities p s , 

V s , 5, and M for various stations along the porous wall can be calcu- 
lated by use of the given set of tunnel stagnation conditions and the 
isentropic flow tables (ref. 6). 

To determine the distribution of the wall-open ratio, r, required 
to provide these desired flow conditions equations (l) and (3) are again 
available; that is, 

m n+ = P a V q sin(S + 7) + cos(5 + 7) + — — cos 7 

m s dx p w p dx sw 

and 

r = mnt 

p jV 

The values of p s , V s , and 5 to be used in equation (l) are those 

determined from the Mach net at the position of the wall in conjunction 
with the stagnation conditions, H q1 and T ol . With the use of the above 

quantities and the given constants of the tunnel, an approximate solution 
of equation (l) can be determined by use of the method of iteration indi- 
cated in the previous section "Estimation of Boundary-Layer Growth," for 
finding values of d&*/dxp with assumed values of wall-open ratio. Of 

course, a solution of equation (l) must be made for many points in the 
flow -generation region in order adequately to define the mass -outflow 
rate along the wall. 

In equation ( 3 ) the values of pj]_ and Vjq are related to the 
plenum-chamber pressure Pj (the value necessary to provide a uniform 

test section) as in step 13 of the preceding section. Also, in this case 
Hjl = H 0 i and T 0 j = T Q p. The quantity o for use in equation (3) is 

determined as mentioned earlier by use of curves of a plotted against 
V s /Vj as functions of M for the material. Approximate values of r 

can now be calculated from the approximate mass-outflow quantities m n t . 
The values of r obtained at this point are only approximate values 
because only assumed or guessed values of r were considered in the 
solution of m n ^.. Thus, an iteration involving r as well as V n t/^s 
must be used as in the preceding section to determine m^ . The final 
value of r is determined by using the approximate values obtained from 
succeeding calculations in performing the next, until the degree of 
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convergence desired is reached. This procedure is used to determine the 
distribution of r for the flow-generation section of the tunnel. 

Since the pressure Pj in the plenum chamber outside the porous 

walls will normally be constant over the length of the nozzle and test 
section and since the stagnation conditions H 0 i and T 0 i for the flow 
in the nozzle are constant, pj and Vj are constant (and equal to Pji, 
Vji) for any particular nozzle design. 

In the -foregoing method, it was assumed that the open ratio of the 
test section was determined by the method presented in the section 
entitled "Selection of the Test-Section Wall Open Ratio." Insofar as 
the design of the flow generation is concerned, any value of r (within 
'reasonable limits) could be chosen for the test section. A different 
choice in the test section r would simply require a different plenum 
chamber Pj necessary to provide uniform flow in the test section. From 
this pressure Pj new values of Pj and Vj, for use in equation (3), 
would be determined. 

Although the flow-generation design method is rigorous in principle, 
it is possible to make serious errors in a tunnel design if the boundary- 
layer growth rates are not known with reasonable accuracy. A nozzle flow 
calculation which will be presented in a later section may be used to 
examine the accuracy to which these quantities must be known. 


APPLICATION OF THE DESIGN METHOD 


In order to evaluate experimentally the effectiveness of the previ- 
ously discussed design method, an approximately 3- by 3-inch two-dimensional 
supersonic porous-wall wind tunnel was designed, constructed, and tested. 

The results determined from using the method for selecting the test-section 
wall-open ratio will be given first and will be followed by a description 
of how the "Flow-Generation Design Method" was used in the design of the 
nozzle . 

It was pointed out in the section entitled "General Design Method" 
that certain design constants or quantities must be chosen; those chosen 
for the design of this approximately 3- by 3-inch wind tunnel are as 
follows: 


Mq 1.278 

T 0 , °F 200 

Hq, ib/sq ft 2086 . 51 + 


(= 29.50 in. Hg) 
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Initial family, v a , supersonic flow expansion rate, deg/in. ... 1.86 

Tunnel height, in 2.954 

Tunnel width, in 5.018 

7, deg 0 

Side walls Solid and parallel 

Top and bottom porous walls Flat and parallel 


The reasons for choosing these particular constants were (l) so 
that the supersonic nozzle and test section resulting from this partic- 
ular design could be tested in an available test apparatus and (2) so 
that the tunnel Mach number and flow angles resulting from proposed model 
tests would be in a range in which other data were available for compari- 
son. Because the boundary- layer growth-rate curves of figure 2 had not 
been determined at the time this tunnel design was made and because no 
other thorough system was known for estimating the wind-tunnel boundary- 
layer growth, a departure from the general design method was made to 
obtain a quantity for m^j to represent the boundary-layer growth mass- 
flow rate . 


The quantity m^j. = 0.008 was chosen to represent the mass-outflow 
rate required in the test section for both the glass side walls and the 
perforated wall. The value of m = 0.008 was derived partially 
through the use of some experimental data of m^i, from tests made in a 

tunnel having 18-percent-open perforated walls, and partially by a guess 
as to how the boundary- layer might grow over a wall with an increasing 
value of open ratio. 


Calculation of Test-Section Wall-Open Ratio 

With the exception of the determination of the boundary-layer growth 
mass-flow rate the design method given in the section on "Method for 
Selecting the Test-Section Wall-Open Ratio" was used for calculating the 
open ratio for the test section of the approximately 3-inch-square tunnel. 

The quantities M 2 , H 2 ^H ol , p q2 , H 2 , p o2 , p 2 , and V 2 for the 
3-inch-square tunnel were determined for the range of 5 2 's given in 
table I. The values found for M 2 , H2/ H ol-> & 2 ’ and v 2 have been 
listed in table I. 

The quantities Pj 2 and Vj 2 used in equation (3) were found by 
first using an approximate value of Pj^ = P3.J to obtain approximate 
values of pj 2 and Vj 2 for computing approximate values of r for 

the range of 62' s in table I. From these data a value of r = 0.47 
appeared best to represent values of 62 between 2° and 3°. 



NACA RM L55J13a 


17 


Next the range of Pj/Hj values in table II were assumed. The 
calculations in table II for Pj, Pjl> Vji, and m n ^. were then made 
with use of the value of r = 0.47* Examination of the data of table II 

P 1 

showed m n ^ = 0.00799 when = 0 . 369 * Since this quantity of 0.00799 

H ji 

was approximately equal to m n ^j = 0.008 the P^ for the tunnel was 
769.87 pounds per square foot. The values of pj 2 and Vj2 for the 
range of &2 were next computed. The curves of figure 4 were used to 
determine values of a. 

Since a constant value of 1 %^ = 0.008 was assumed for all values 
of 62 no iterative solution was required to determine 111^2 from equa- 
tion (l). Table I shows the values determined for The values of 

V 2 /Vj , a, and r (shown in table I) were then found for the range of 
6 ‘ 2 f s chosen. The results of these calculations showed that a value of 
r = 0.47 should be a reasonable average for use in the test section. 


Since the method described for selecting the test-section wall open 
ratio indicates that the boundary layer should be estimated by the use 
of curves similar to those of figure 2 , calculations were made to find 
what values of wall-open ratio would result when using this system. 

Table III shows some of the quantities used in making these calculations. 
Notice that the values of the open ratio given in table III are not appre- 
ciably different from those of table I for values of 62 = 2 . 57 ° or above. 
This approximate duplication occurs because the boundary-layer growth 
estimates used for the two cases form only a small percentage of the 
total wall outflow. The greatest influence of the boundary-layer growth 
occurs at the low values of bg. 


The open-ratio values given in tables I and III for 62 = 0 can be 
used as an index to show what values of open ratio were used in the itera- 
tion solution to find the other values of r in the table. Remember the 
value of r assumed to find Pj (62 = 0) should, in the final solution, 

be approximately equal to the r selected for the test section. 

Although the total pressure losses across an oblique shock wave are 
included in the calculation of tables I and III, no appreciable errors 
would have been caused if these shock losses were neglected. If the data 
of table III had been used in selecting the test-section open ratio, a 
value of r slightly lower than 0.47 would have been chosen. A sample 
calculation is given in appendix A to show how the boundary layer was 
estimated for table III. 
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Design of Nozzle 

According to the method as stated in the section on flow generation, 
a supersonic characteristic flow net was constructed for a 1.86° per inch 
first family v a expansion rate (ref. 9)- The characteristic net was 
completely contained between parallel boundaries h/2 apart, and the 
characteristic lines were initiated at the boundary (chosen to represent 

1 O/TO 

the porous wall) for each ~ L '° D v a family expansion. All characteris- 
tics intersecting the wall boundary from the family were assumed 

to have no downstream wall reflection, and thus, were terminated at the 
wall intersection point. This type of construction will result in the 
development of a network which will yield a parallel flow in the minimum 
distance. With the nozzle characteristic network complete, it was possi- 
ble to determine the flow expansion angles along the tunnel walls. 

Table x IV shows the expansion angles determined for the design of the 
approximately 3-inch-square wind tunnel. These expansion angles were 
used to compute 6 S = (v a - v^), M s , p s , and V s , at each chosen wall 
station for the given tunnel stagnation conditions. 

It was pointed out in the "Flow Generation" section, that equa- 
tions (l) and (3) should be used to calculate the required open ratio r 
of the wall, for wall-flow conditions, as determined for the character- 
istic network, and for the values of Pj-^ and as determined from 

the conditions of the empty test section. This procedure was followed 
in the design of this nozzle with the exception that the boundary-layer 
growth terms were evaluated differently. The reason for making the 
change in estimating the boundary layer has been explained earlier. The 
boundary-layer mass-flow rate distribution term, m^^, in this design 

was estimated by first neglecting boundary-layer growth in the flow gen- 
eration region and computing an r distribution with the use of equa- 
tions (l) and (3) • The values of Pj-^ and Vjj_ used in the calculation 

were the values determined in the previous section "Selection of the Test- 
Section Open Ratio" using a test section m^^ = 0.008 , 6p = 0, and 
r = 0.47. 

Examination of the resulting r distribution showed that the tunnel- 
open ratio increased rapidly with increasing downstream distance until 
the nozzle station 1.5 inches downstream of the entrance of the flow- 
generation region was reached. At this station, the nozzle open ratio 
was approximately that of the test section; from this point on, the open- 
ratio distribution became nearly constant throughout the remainder of 
the flow-generation region until the last few stations were reached. 

These data led to the boundary- layer mass-flow estimates m n ^2, given 

in table IV. The boundary-layer growth mass flow m^j was assumed to 
increase linearly with distance (from 0 to 0.008) until the 1.5-inch 



NACA RM L55J13& 


19 


station was reached, and then to remain constant for the following down- 
stream stations. This approximation of the boundary-layer growth esti- 
mate was added to the outflow necessary for supersonic expansion to give 
the total outflow rate m^ at the various tunnel stations. These data 
and the porosity characteristics a taken from data of figure 4 were 
inserted into equation (3) to compute the required open-ratio distribu- 
tion of the tunnel. The results of the calculations are tabulated in 
table IV. A sample calculation for this table is given in appendix B. 
The flow-generation region of the nozzle ends at station 3 *79* The open 
ratio for all stations downstream of this point must be the preselected 
test-section open- ratio value, 0.47* 


Construction and Tests 

The wall-open-ratio distribution for the wind tunnel was achieved 
by drilling two 0. 060-inch-thick aluminum-alloy plates with holes ranging 
in diameter from 0.016 inch to 0.046 inch. Figure 5 shows a photograph 
of one of the walls. The holes were drilled in rows across the tunnel 
width, with all holes in each row having the same diameter. The rows 
were spaced approximately 0.030 inch apart. Each row was drilled so that 
the total open area required for a 0.030 by 3.018 strip would be given by 
the holes. Although considerable care was taken in the drilling of the 
walls, the final wall-open-rat iO distribution as determined by measure- 
ments was slightly different, in general, from that computed in the 
design. The data of figure 6 show a comparison between the measured and 
design open-ratio values. In order to show what these construction 
errors mean in percent of open ratio, the data of figure 7 are given. 

When the measured values of open ratio (fig. 6) were used to com- 
pute the flow field of the nozzle, the data of figure 8 were obtained. 
(The system used in calculating the nozzle flow, with the wall-open ratio 
given, will be explained later.) If the data of figure 7 are used in 
conjunction with the data of figure 8, it is possible to note at some 
points the effect of errors in wall construction on the flow generation. 
For instance, the 11-percent error in wall-open ratio at station 0-51 
can be seen to cause a rapid change in Mach number occurring at sta- 
tion 0.51 of figure 8(a). This error introduced at station 0.51 also 
affects the flow at the center-line station 1.2 and the wall station 2.1 
as shown in figures 8(b) and 8(a), respectively. The next downstream 
center-line station which could have been affected by the disturbance 
resulting from this upstream open-ratio error is station 3*21. .Notice, 
however, that no appreciable flow disturbance shows up at this point. 

This result is believed to occur because the turning angle provided by 
the wall at station 2.1 on account of the presence of the disturbance 
itself is of the correct .magnitude and direction to cancel the incoming 
disturbance. Since the downstream flow is influenced by errors in open 
ratio made upstream, it is difficult to establish from the data of fig- 
ure 8, the exact effect of an open-ratio error at a given station. 
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Nevertheless, it can be stated that the open ratios as measured for the 
porous wall are sufficiently small to provide an approximately uniform 
center-line flow downstream of station 2.6. 


When the tunnel was put into operation, center-line pressure meas- 
urements were taken to determine the degree of flow uniformity. The data 
of figure 9 show the Mach number variations through the nozzle and test 
section for three different settings of plenum-chamber pressure. On this 
figure, the circles represent the experimental data for the design condi- 
tion (M = 1.278), and the solid line indicates the designed Mach number 
distribution i The experimental data show that the flow expansion in the 
nozzle started at least 0.2 inch upstream of station zero. This upstream 
expansion, which is quite common in both perforated- and slotted-wall 
nozzles, is believed to be due to a decrease in the boundary-layer thick- 
ness as the perforated portion of the nozzle is approached. 


The general overexpansion occurring at the center line in the down- 
stream portion of the flow-generation region (stations 1 to 2.6) is prob- 
ably caused by the use of an inaccurate boundary- layer growth rate in the 
design. The assumption made in design that m n bZ i* 1 the flow-generation 
region increases linearly from 0 at station zero to 0.008 at station 1.5 
was probably considerably in error because (l) figure 2 shows that d5*^dXp 

decreases as V n t/V s increases; thus the actual P S V S term would 

' dX-r, 


decrease at first with increasing station in the nozzle and be quite small 
for the high wall outflow rates; and (2) the boundary-layer growth on the 
side walls is approximately constant and becomes the governing component 


of “nbl 


when P S V S is small. 

dx p 


The walls were designed to be flat between the wall stations zero 
and II.5. Surface measurements, however, revealed that wall waves of 
±0.003 inch existed on each wall. The exact effects of these surface 
contours on the flow could not be definitely determined; however, these 
wall variations must have influenced the flow to some degree. 

The experimental data taken at Mach numbers 1.1 and 1.2 (fig. 9) 
are included to show the type of flow distribution which was obtained 
at off-design Mach numbers. A comparison of these data with those taken 
at the design Mach number shows that much more severe flow variation 
occurred at the off-design Mach numbers. 

Insofar as the testing of models is concerned, the flow variations 
in this 3- by 3-inch tunnel are sufficiently small to provide a suitable 
test region. For example, at the design Mach number (M = 1.278), the 
maximum Mach number variation along the center line between the sta- 
tions 2.5 and 12.0 is AM = ±0.007 - The rapid flow generation, which 
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was achieved in less than one tunnel height, and the extended run of 
fairly uniform flow, are most encouraging in view of some of the uncer- 
tainties involved in estimating the boundary-layer growth, and of the 
physical precision of the fabricated wall. 


Calculation of a Wind-Tunnel Flow Field 

As an outgrowth of some of the principles used in the nozzle design 
described in this paper, an effective method has been employed for cal- 
culating the wind-tunnel flow field. This method primarily involves the 
finding of the flow conditions at points along the wall and using these 
flow conditions to determine the wind-tunnel flow field. The discussion 
given here will cover only the method for finding the flow conditions at 
one station along a porous wall. Procedures are already established for 
determining the stream flow once the flow at the boundary is known 
(ref. 9 ). 

To determine the values of M and 5 at a given point on the wall, 
it is necessary to know the flow conditions just upstream of the point 
to be computed. This requirement will in general, make it necessary to 
start the calculations at the tunnel minimum, where the flow is uniform 
and M = 1.0. It is also necessary to specify the wall boundary- layer 
growth, the wall convergence, wall -porosity characteristics a and open 
ratio, and plenum-chamber pressure. With these quantities given and with 
the aid of equations (l) and ( 3 ), the flow removal rate through the porous 
wall at any station providing for supersonic expansion in a tunnel can be 
found by solving the following expression: 


m ex Pjl\jl ra “ m nbZ - m nt “ m nbZ - Ps^s sin(5 + 7) ( 5 ) 


For the special case of parallel walls, equation ( 5 ) becomes: 

“^x = Pjl^jl ra “ m nbZ = Ps v s s ^ n S ( 6 ) 

To determine the Mach number resulting from this outflow, first find 
the flow angle 6 p t at the point. This flow angle may be used in con- 
junction with the upstream flow conditions of 5 a ^ and v ah to find the 

expansion angle at the point. The following equation shows how these 
quantities are related: 


v pt = v ah + ( 5 pt " 8 ah) 


( 7 ) 


The Mach number can be determined from the expansion angle v 1 
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An examination of the variables in equation (5), rewritten in the 


form 


P S V S sin(6 + 7 ) = -PjlVjiro - P S V S ^ cos(S + 7) + ^ cos 7 (8) 


dx r 


w p dx sw 


shows that to determine the value of 5 at a point from this equation, 
it is necessary to know the Mach number at the point because the quan- 
tities p s V s , d6*|dXp, and a are functions of M. Since M is one 

of the wall quantities to be found from the final solution, equations ( 5 ) 
and (7) must be solved by an iteration process. If, in solving for the 
outflow or 6pt, the stream quantities p, V, M, d6*|dXp, and a just 

upstream of this point are used for the values, of p, V, M, d5*/dx p , 
and a at the point, these quantities will serve as a good first approx- 
imation for the iteration. The Mach number determined from this first 
calculation of v p t may now be used for making the second approximation, 

and the results of the second in computing the third, etc. A more rapid 
convergence for the iteration will occur, however, if the average value 
of Mach number resulting from the two next preceding approximate calcu- 
lations is used for recomputing Mp t and 6 pt . The most important use 

for this wind-tunnel flow-field calculating system is that it makes pos- 
sible calculations of the stream flow conditions at any off-design Mach 
number, and thus should be useful to one who wishes to design a transonic 
tunnel to operate over a range of Mach numbers. 

The circles shown in figure 8 are calculated for the 3- by 3-inch 
tunnel by use of the design specification and the measured open-ratio 
distribution. These data were determined by calculating the flow field 
produced by each lateral row of holes. The average row spacing was 
0.03015 inch. Normally, only three iterations were required to satisfy 
the conditions at the wall to a Mach number accuracy of ±0.00004. The 
Bell Telephone Laboratory X-66744 relay computing machines were used to 
make these calculations. For the reader's convenience, a- sample calcu- 
lation of the data of figure 8(a) has been included in appendix C. 


Discussion of the Wave-Reflection Problem 

Since a design method has been given for calculating the open ratio 
needed in a porous -walled test section to minimize shock-wave reflection 
ability of the wall, some discussion is in order to indicate the phenom- 
ena which can be encountered during the practical use of such walls. 

In the design method, it was assumed that-- the boundary- layer growth 
rate in front of the shock wave would be constant (in test-section region) 
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It was also assumed that the boundary-layer growth rate would be con- 
stant for a given value of 62 behind the oblique shock wave. Although 
these assumptions were necessary to simplify the design problem, it was 
known that these conditions did not exist in practical wind tunnels. For 
instance, the boundary- layer thickness over a porous wall was found in 
reference 3 to change rapidly in the vicinity of the oblique shock wave. 
The variation in boundary- layer thickness upstream (from shock) resulted 
because the increased pressure behind the shock traveled upstream through 
the boundary layer. (Ref. 10 shows this phenomenon occurring over solid 
walls.) The effect of this pressure behind the shock causes the upstream 
boundary layer to grow and reach a peak thickness just ahead of the model 
shock-intersection point and results in the generation of a compression 
region which in some cases generates a shock in the stream. Immediately 
downstream of this point, the boundary layer thins because of the pres- 
sure rise occurring across the model shock. The thinning of the boundary 
layer increases the effective area ratio of the tunnel and thus creates 
an expansion wave. The strength of the expansion is magnified to some 
extent at this point by the porous walls because of the ability of the 
wall to remove flow from the tunnel. 


It is believed that if the boundary layer over the walls of the 
tunnel is made thin, as is possible by designing the tunnel to have a 

the effect of the boundary 

U-X-n 


value of 7 much greater than arc tan 


layer on the wave reflection of the model shock would be much less severe 
and would induce only minor flow variation in the stream. 


It should also be recognized that, if the boundary- layer growth 
rates were known in the vicinity of the point where a. given strength 
model shock wave intersects the wall, equations (l) and ( 3 ) could be 
used, theoretically, to calculate a distribution of open ratio along 
the wall in vicinity of the shock which would eliminate the reflection 
of the shock. Since this solution for the wave-reflection problem 
accounts for .only one flow condition in the tunnel it appears to be 
impractical for use in wind-tunnel testing. More research will be 
required before practical solution to the wave-reflection problem can 
be demonstrated. 


CONCLUSIONS 


The following conclusions summarize the wind-tunnel-design method 
and experimental, information described in this report: 

1. An effective method is presented for the design of transonic 
porous-wall wind tunnels which at present is known to be capable of 
providing a flow suitable- for testing at the design Mach number. 



2k 


NACA RM L55J13a 


2. The design method permits a supersonic tunnel to be designed to 
have, within limits, any arbitrarily selected rate of flow expansion in 
the nozzle or any desired center-line Mach number distribution. 

3. Flow generation for a design Mach number of 1.278 can be achieved 
in. less than one tunnel height with test-section Mach number variations 
within ±0.007. 

4. A method is given for calculating the open ratio of the wall 
needed to maintain a constant pressure field behind an oblique shock 
wave . 


5. The constant value of open ratio in the test section, even though 
chosen by the design method, is not expected to prevent shock-wave- 
reflection disturbances in the presence of any appreciable wall boundary 
layer. 

6. Model shock-wave reflections from walls having thin boundary 
layers are expected to be quite weak when the open ratio of the test 
section has been chosen through the use of the design method. 

7. The quality of the design depends to a considerable degree on 
how accurately the boundary-layer growth rates over the wall and the 
porosity characteristics of the porous wall have been determined. 

8. The wind-tunnel flow field in a porous-wall tunnel can be deter- 
mined for any transonic Mach number by the specific calculation procedure 
given in this report . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 6, 1955* 
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APPENDIX. A 

BOUNDARY-LAYER ESTIMATION 

In order to illustrate how the boundary- layer growth was estimated 
in table III, the following calculations were made for 62 = 0 (the uni- 
form flow condition in the test section): 

Constants given for the design of the 3- by 3-inch tunnel are as 
follows : 

Mi = 1.278 
H ol = 2,086.54 
T 0 i = 200° F 
p 0 i = 0.0018434 

7 = 0° 
a©! = l,26l 

p-L = 0.000909 

Vx = 1,399.2 

Boundary- layer growth-rate curves, figure 2 
Porosity-characteristic curves, figure 4 
Height of solid side walls, 2.95^ in* 

Width of perforated walls (top and bottom), 3*018 in. 

Boundary- layer growth of side walls, - = 0.002 

ax sw 

The following assumptions were made: 

1. Assume that the boundary- layer-growth curve for M = 1.22 in fig- 
ure 2 applies for the perforated wall at M = 1.278. 
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2. Assume that the required test-section open ratio will be 0.47. 
(The results of table I indicated that this value should serve as a good 
first approximation.) 

3. Assume as in the section on "Estimation of Boundary- Layer Growth" 
that the quantity 0.47/0. 4l times the boundary-layer growth rate for 
4l-percent open material gives the boundary-layer growth rate for a 
47-percent open wall. 

To find the boundary- layer growth mass flow, we must find a solu- 
tion for equation (4) 


m nt2 = P2 V 2 sin ( 5 2 + 7) + P2 V 2 ^ - 5 * C °^ 


dS* cos(52 + 7) , h-sw dS* cos 7 


W r 




S W 


For parallel walls and free-stream conditions in the test section 
(6 = 0°), P 2 V 2 = PiVp. Therefore, 


%itl - PlVi^sin 0 ° + ^ ^ 


In order to evaluate the quantity 
To start the solution, first assume 


(mnti)i = 0.000909 X 1399.2^0 + 0 


d5*^dXp, an iteration must be used. 

= 0 and 4^- = 0.002 then 
dx p dx sw 


+ 1954 q 

3.018 


.002J = 0. 


00249 


Subscripts numbers outside parenthesis denote the number of the iteration 
step. 



(mntl)l _ 0.00249 
p l 0.000909 


2.74 


(Vntl 4 = -Si2L = 0.00196 
Vl 1399.2 


This is the outflow- velocity ratio required to compensate for 
boundary -layer growth on the side walls. 
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Now by using 
first approximate 


( v ntl)l 

Vl 

value of 


= 0.00196 and the curves of figure 2, 
d6*/dxp will be found to be 


the 


d5* 

dx p 


0. 00625 x 


0. 007165 


Pl v l 


d5* 

dx p 


0.009113 


With this term evaluated, the second approximation to the total mass- 
outflow quantity is 


( m ntl)2 = (0 + 0.009113 + 0.00249) = 0.11603 


and the velocity ratio 


( v ntl)2 _ 12.76 
V x " 1399.2 


0.00912 


Since (Vntl^/v'l does not equal (^ntl)l^l an °ther choice of 
outflow velocity must be tried- For rapid convergence, one-half the 

previous value of P]_Vq or is applied for a third 

P 

approximation 


( ra ntl)3 = 0.004557 + 0.00249 = 0.007047 

. fa -Us . i-r» = 0.00555 
Vi 1399.2 


— = 0.00292 X = 0.00335 
dx p 7 0.41 


for 


( V ntl) 3 / V l’ 
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Pl v l T~ = °* oo42 57 

dXp 


( m ntl) 4 = °*0°^ 2 57 + 0.00249 = 0.006747 
( v ntl)4 7-42 


1399.2 


= 0.00531 


for 




= 0.00310 x 111 = 0.00355 

dXp 0.41 


W = - 1 * 1 . . 0.0051.7 

V x 1399.2 


The desired value of 

^ = 0.00355 and 0.00335. 
dx p 

a more rapid convergence 


d&* j dXp must 
By using the 


fall between the values of 
average of these two values for 


d&* = 0.00355 * 0.00335 = 0 00 , 45 

dXp- 2 

= 7.56 _ _ 0.0054 
Vi 1599.2 


for ( V ntl)6/ V l> 

^ = 0.00301 x 2iiil = 0.00345 

^ 0.41 

Since the value of = 0.00345 is the same value as that used 

dx p 

in the preceding step, this quantity 0.00345 slug/ft 2 -sec represents the 
boundary-layer growth rate over the perforated walls in the test section. 
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APPENDIX B 

DETERMINATION OF OPEN RATIO IN NOZZLE 


A sample calculation has been prepared to show how to find the 
open ratio for a station along the wall in the flow-generation region 
(table IV), By using the constants listed in section "Application of 
the Design Method" and proceeding as indicated below, the quantities in 
table IV for the station 1.2 were found. 

The values of v a and in table IV were read from curves which 

were determined from the nozzle characteristic flow network. 

At the wall station 1.2, 


v a = 2.233° and = 0 . 269 ° 
v s = 2.233 + 0.269 = 2 . 502 ° 

6 S = 2.233 - 0.269 = 1 . 964 ° 

From the isentropic-flow tables, the value of M for v = 2.502 is 
M = 1.1553* For this Mach number, the following stream quantities were 
determined by use of the constants H 0 p = 29-50 in. Hg(= 2,086.5 lb/ft^), 
and T 0 i = 200° F: 


p ol 


Hpi _ 2086.5 

RT 0l 1715 x (200 + 460) 


0.0018434 


a ol = ^9*ljT^ = 49.1^200 + 460 = 1,261 ft/sec 

a s = 1,120.3 ft /sec for M = 1.1553 

V s = M s a s = 1.1553 X 1120.3 = 1,294.3 ft /sec 


The quantities pjj. and Vjp to be used in the following equation, 
equation ( 3 ) of the text, 


m nt 

pj v j° 


r 
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are the values determined for 62 = 0 in table I. These values are 
pjl = 0.000904 and Vjl = 1403-8. Now having determined the value 
of Vj-|_, the quantity V s / Vj for the station 1.2 of table IV' is found 
to be • x 


_ 1 gg4 . i | = 0.9220 
Vji 1405.8 


By use of the curves of figure 4, the value of a corresponding to 
the stream Mach number (M = 1.1553) and velocity ratio V s jVj = 0.9220 
was found to be 0.1370. 

The total wall mass -outflow rate m n -t is obtained by the use of 

equation (1) with the exception of the terms including the boundary-layer 

growth. The mass-outflow rate 1%^ necessary to compensate for the 

boundary- layer growth was found by the system stated in the section 

entitled "Application of the Design Method," where 1%^ was assumed to 

vary linearly from a value of zero at station zero to a value of 0.008 

at station 1.5. Thus, the quantity m n ^x = °- 0Q Q ^ lr2 = 0.0064. 

-L • 5 

Now to obtain from equation (l) the outflow component providing for 
supersonic expansion 

mgx = P s V s (sin 6) for parallel walls 


Hi- = 0.6510 for M s = 1.1553 
p ol 


p s 



= 0.6510 X 0.0018434 = 0.00102 


Thus, the value of the supersonic expansion mass outflow is 


mex = 0.00102 X 1120.3 sin 1.964° = 0.04524 


and 


m n -t = n^nbl + = 0.0064 + 0.04-524 = 0.05162 
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therefore. 


r 


m nt 

p ji v ji° 


0.0^162 

0.000904 x 1403.8 x O.1370 


0.2968 or 29.68 percent 


This value of r = 29-68 percent is the open ratio the wall must 
have to satisfy the desired flow conditions at the 1.2-inch station. 
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APPENDIX C 

DETERMIMTION OF LOCAL MACH NUMBER AND FLOW 
ANGLE AT A POINT ON THE WALL 


In order to find the Mach number and flow angle at a given wall sta- 
tion in the nozzle of the 3- by 3-inch tunnel constructed by the presented 
design method, the following sample calculation is made: 

Given quant it ie s : 

M = 1.0 upstream of point to be calculated 
v s = 0° upstream of point to be calculated 
5 S = 0° upstream of point to be calculated 
P ol = 0.0018434 

7=0° 

T ol = 200° F 

Station = 0.03015 in. (station for which this calculation is made) 

m = 0.00016 slug/ft 2 -sec (estimated by system described in the 
section "Design of Nozzle" ) 

Vj]_ = 1403.8 (from table IV) 

Pj3_ = 0.000904 (from table IV) 

Porosity-characteristic curves of figure 4 

Wall-open ratio curves of figure 6 (r at station 0.03015 in. is 0.0047) 

The Mach number at the given station, 0.03015, may be determined by 
finding the solution of equations (5) and (7) at this point. Equation (5) 
is given as 

“ex = (Pji v jl ra " m nbz) = PsVs sin ( 5 + 7) 
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for the case where 7=0 equation (5) may be rewritten as 

“ex = Pjl v jl ra " “nbi = P s v s si * & 
and equation (7) is given as 

v pt = v ah + ( 5 pt " 6 ah) 

The unknown quantities which will be determined are a, p s , V s , and 6. 

As was explained earlier a, p s , and V s are all functions of the 
Mach number or v pt> the quantity which is to be determined; thus, the 
solution must be found by the use of an iteration process. 

To do this, first assume that the velocity at the station being cal- 
culated will be the same as that at the upstream station and find the 
Mach number resulting from this solution. By the use of this upstream 
point, 


81 

v s = M s — a o = i:l 5Q-3 (for M = 1.00) 
a 0 


Vs_ 

v jl 


1130-3 = 
1403.8 

a = 0.279 


0.820 


Then, 


rn ex = 0.000904 X 1403.8 X 0.0047 X 0.279 - 0.000160 = 0.001504 


P S 


Ps 

p ol 


p ol = 0.6339 X 0.0018434 = 0.001169 (for 


M = 1.000) 


6 = 6 . 


'pti 


= arc sin 


m ex 

P V 
H s s 


= arc sin 


0.001504 


O.OOH69 x 1150.5 


= 0.06409° 
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The subscript pt}_ is being used to denote the first partial results at 

the point for which the solution is sought. With increasing iteration 
steps , the partial results will be denoted by pt 2 , ptj, etc. 

v ptl = v ah + (& pt - &ah) = 0 + (O. 064 l ' °) = °- 06 ^ 10 


The Mach number resulting from v ptl (M = 1 . 015 ) is used to make 
the following second approximation: 

for Mpt x : 


V s = Mp t ap tl = 1.015 X 0.9H0 X 1261 = 1165.7 

Vs_ _ = o.8290 

Vjp 1405-8 

a = 0.2675 

me X = (0.000904 x 1405.8 X 0.0047 x 0.2675) - 0.00016 = 0.001455 
for Mpt pi 


P s = 0.001155 

6^p = arc sin 0 - 001 435 = 0.06106° 

p 2 1165.7 x 0.001155 

v pt2 = 0 + (O.6106 - 0) = 0.06106° 

Mpt 2 = 1.01250 

Similarly the third approximation: 
for M pt 2 : 


V s = M p t 2 apt 2 = 1.0125 X 0.9111 X 1261 = 1165.5 
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= 0.8286 

VJ1 1403.8 
a = 0.268 

nigx = (0.000904 x 1403.8 x 0.0047 x 0.268) - 0.000160 = 0.001438 


for M pt 2 : 


P S - Pol - 0.001156 


0.001438 ' o 

6 P t, = arc sin = 0.06112° 

5 1163.3 X 0.001156 

v p t 5 = 0 + (0.06112 - 0) = 0.06112° 

Mp tj = 1.01252 

The value of Mp-tjJ 6 ives the Mach nutnher at station 0.03015 to an 
accuracy of ±0.00002. 
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TABLE II.- TABULATION OF SOME VALUES USED TO FIND THE PLENUM-CHAMBER PRESSURE FOR A 
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TABLE III.- SOME QUANTITIES USED IN CALCULATING THE REQUIRED OPEN RATIO FOR THE 3-INCH 






TABLE IV.- SOME QUANTITIES USED IN CALCULATING THE OPEN-RATIO DISTRIBUTION IN THE 
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Outflow velocity ratio, 
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Figure 2.- Curves showing how d5*/ dx p varies with the outflow velocity 

ratio for several Mach numbers over a 4l-percent-open perforated 
specimen. 
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Outflow velocity, V n ^ , ft/sec 

Figure 3.- Curves describing how the friction coefficient varies with 
outflow velocity for a 4l -percent-open perforated specimen. 
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Figure 4.- Discharge coefficients for a *4-1 -percent-open perforated speci- 
men having an incoming boundary -layer displacement thickness of approxi- 
mately 0.002 inch. 
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